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PREFACE 

Within  the  framework  of  the  NBS  Institute  for  Materials 
Research  the  area  of  standard  reference  materials  Is  a  broad 
and  Important  one.  Including  the  preparation,  characterization 
and  distribution  of  a  wide  variety  of  materials  In  such  diverse 
fields  as  metallurgy,  polymers  and  Inorganic  materials.  In 
carrying  out  such  a  program  there  Is  much  Interaction  with 
representatives  of  Industry  and  science,  beginning  with  dis- 
cussions as  to  which  primary  standard  materials  will  do  most 
to  advance  technology,  the  furnishing  of  materials  and  fabri- 
cation of  samples,  and  the  characterization  and  certification 
of  the  materials  by  cooperative  efforts.     The  many  groups  par- 
ticipating In  a  standards  program  are  very  Interested  In  de- 
tailed Information  on  specific  aspects  of  the  program  —  but 
to  date  there  has  been  no  publication  outlet  for  such  written 
discussions . 

To  meet  this  need,  NBS  Miscellaneous  Publication  260  has 
been  reserved  for  a  series  of  papers  In  the  general  area  of 
"standard  reference  materials".    This  series  will  present  the 
results  of  studies  and  Investigations  undertaken  within  the 
Institute  for  Materials  Research  with  emphasis  on  the  prepara- 
tion and  characterization  of  standard  reference  materials. 
This  subject-oriented  series  will  provide  a  means  for  rapid 
dissemination  of  this  detailed  Information  and  we  hope  will 
stimulate  the  use  of  standard  reference  materials  In  science 
and  Industry, 

W.  Wayne  Melnke,  Chief 

Office  of  Standard  Reference  Materials 
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MOSSBAUER  SPECTROSCOPY  STANDARD 
FOR  CHEMICAL  SHIFT  OF  IRON  COMPOUNDS 


Jon  J.  Spijkerman,  David  K.  Snediker^ 
Flllmer  C.  Ruegg^  and  James  R.  DeVoe 

Institute  for  Materials  Research 
National  Bureau  of  Standards 

ABSTRACT 

The  preparation,  calibration,   and  use  of  the  Standard 
Reference  Material  for  chemical  shift  of  iron  compounds  in 
Mossbauer  Spectroscopy  is  described.     This  standard  is  a  pro- 
perly oriented  single  crystal  of  sodium  pent acyanonitrosyl- 
ferrate(ll)  dihydrate  ( sodium  nltroprusside ) .     Primary  stan- 
dards were  calibrated  with  a  high  accuracy  optical  Mossbauer 
spectrometer  and  secondary  calibrations  were  made  with  a  tan- 
dem Mossbauer  spectrometer.     The  midpoint  between  the  two 
absorption  peaks  in  the  spectrum  at  25.0  "^C  provide  a  useful 
velocity  scale  calibration  of  the  spectrometer.  Descriptions 
of  the  spectrometers  used  and  the  error  analysis  associated 
with  the  data  are  given.     Suggested  format  for  reporting 
Mossbauer  spectra  and  their  parameters  is  also  presented. 


Key  words:     Mossbauer,   Spectroscopy,   standard,  differential 
chemical  shift,   quadrupole  splitting,  random  error,  systematic 
error,  sodium  nltroprusside,  iron  compounds,  optical  Mossbauer 
spectrometer,  tandem  Mossbauer  spectrometer,  velocity  calibra- 
tion. 
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1.  INTRODUCTION 

While  the  Mossbauer  effect  was  first  used  to  resolve  a 
number  of  problems  in  physics  it  has  now  found  a  definite 
application  in  chemical  structure  analysis.     The  Importance 
of  Mossbauer  spectroscopy  to  chemical  structural  analysis  was 
recognized  in  I962  but  received  little  attention  until  1964. 
The  early  efforts  made  to  explain  the  parameters  obtained  in 
Mossbauer  spectroscopy  were  qualitative  or  at  most  in  a  semi- 
quantitative sense.     Only  since  I965  have  attempts  been  made 
to  be  quantitative.     Such  attempts  have  brought  about  realiza- 
tion of  a  serious  shortcoming  with  data  published  in  scientific 
periodicals.     The  Mossbauer  chemical  shift  for  a  given  series 
of  compounds  is  rarely  reported  with  respect  to  the  same  source 
(see  figure  l).     Consequently,  normalization  of  data  taken  by 

different  laboratories  is  necessary.     One  additional  diffi- 
culty has  been  found.     Even  if  the  source  used  in  two  differ- 
ent laboratories  is  nominally  the  same,,   serious  descrepancies 
occur  in  the  chemical  shift  reported.     This  presumably  is  a 
consequence  of  the  method  of  preparation  of  the  sources  which 
usually  are  made  by  diffusing  a  radioisotope  into  a  metal 
substrate.     Proper  annealing  is  very  important  and  doing  this 
reproducibly  to  give  the  same  structural  environment  is 
difficult . 

The  most  practical  solution  is  to  measure  the  chemical 
shift  of  an  unknown  compound  with  respect  to  a  reference  com- 
pound.    This  was  anticipated  by  the  attendees  of  a  symposium 
on  Mossbauer  Effect  at  the  1964  Gordon  Research  Conference  on 
inorganic  chemistry.     The  compound  sodium  pentacyanonitrosyl- 
ferrate(ll)  dihydrate  (sodium  nitroprusside )  was  suggested  at 
this  conference.     Other  materials  were  also  suggested.  A 
particularly  interesting  one  appeared  to  be  iron  foil,  but 
studies  at  NBS  have  indicated  that  its  reproducibility  with 
respect  to  both  production  of  the  standard  and  stability  of 
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Figure  1.     Chemical  shift  of  sodium  nltroppusslde  for 
different  source  matrices. 


It  was  necessary  to  develop  a  philosophy  for  the  produc- 
tion and  calibration  of  this  standard  which  would  maintain 
Its  usefulness  during  the  foreseeable  future.     The  system  de- 
cided upon  Is  like  that  of  many  of  the  standards  developed  at 
NBS.     Ten  primary  standards  were  made  from  the  best  obtainable 
single  crystal  of  sodium  nltroprusslde .     These  crystals  were 
then  calibrated  with  a  spectrometer  that  has  very  high  accu- 
racy and  precision.     Fifty  standards  were  then  randomly  sara- 


pled  from  a  lot  of  500.     The  sample  of  fifty  was  then  cali- 
brated against  the  primary  standards  using  the  NBS  tandem 
spectrometer.     The  ten  primary  standards  will  he  measured 
periodically  with  the  high  accuracy  spectrometer  to  detect 
any  possible  changes  in  the  chemical  shift  as  a  function  of 
time.     Any  subsequent  lots  of  standards  will  be  calibrated 
against  the  ten  primary  standards.     In  addition,  if  standards 
for  compounds  other  than  those  containing  iron  should  be  de- 
sired, these  standard  reference  materials  would  be  calibrated 
against  the  primary  standard  for  iron. 

This  report  describes  the  theory  and  procedures  for  cal- 
ibration and  use  of  the  Standard  Reference  Material  for  chem- 
ical shift  of  iron  compounds.     It  has  not  been  found  practi- 
cable to  include  all  of  the  detail  on  the  single  crystal  pre- 
paration, on  the  Mossbauer  spectrometers,  or  on  the  computer 
programs  used  to  resolve  the  spectra.     Pertinent  references 
on  these  subjects  are  giveno 

It  is  important  that  this  new  and  powerful  spectroscopy 
receive  the  benefit  of  standardization  so  that  the  spectrum 
parameters  can  be  used  with  confidence.     In  conjunction  with 
this,  a  format  for  reporting  Mossbauer  spectra  and  for  re- 
porting spectrum  parameters  along  with  their  nomenclature  is 
suggested . 
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2.     PREPARATION  AND  rFSCRIPTION  OF  THE  STANDARD 

A.  Preparation 

Sodium  nitroprusside^  ^R^[Fe{m)^m]'2E^0,   is  prepared 
by  the  action  of  ^0%  (vol)  nitric  acid  on  sodium  ferrocyanide 
or  ferricyanlde   [5],  or  by  the  reaction  of  nitride  ion  on 
either  sodium  ferrocyanide  or  ferricyanlde  in  acid  solution 
[^.5,6]. 

Single  crystals  were  grown  from  a  saturated  aqueous  so- 
lution of  ACS  reagent  grade  material  [7]  at  a  pH  of  6.5  -  6.8^ 
using  glass  equipment.     Thin  platelets  1  cm  x  1  cm  x  .O78  cm 
were  cut  from  the  single  crystals.     The  platelets  were  cut 
with  a  string  saw  along  the  100  plane  within  t  2°  .     The  re- 
quired thickness  of  0.0775  -  O.OO3O  cm  with  a  parallelism 
within  0.001  mm  and  a  surface  finish  of  20  microns  was  achiev- 
ed by  polishing  the  platelets.     A  stock  of  500  platelets  was 
produced  by  the  Isomet  Corporation  for  the  NBS  standard  [8]. 

B.  Chemical  Description 

Sodium  nitroprusside  Is  stable  at  room  temperature  under 
normal  laboratory  conditions  [9].     The  compound  with  two  water 
molecules  of  hydration  is  stable  up  to  100° C.     At  higher  tem- 
peratures (near  450°C),  the  salt  is  transformed  [10]  to  sodium 
ferrous  ferrocanlde,   sodium  ferrocyanide^  nitric  oxide  and 
cyanogen.     At  higher  temperatures^  nitrogen  is  evolved.  To 
determine  the  stability  of  the  standard  reference  crystal  in 
a  laboratory  atmosphere,   a  platelet  was  exposed  for  a  prolong- 
ed period  of  time.     The  relative  humidity  was  50-70^  and  the 
temperature  25  ±  2°C.     No  change  in  weight  of  the  platelet  was 
observed  within  10"'^  g  after  a  six  month  period.     For  stability 
with  respect  to  the  Mossbauer  spectrum,   see  section  3. 

The  solid  salt  is  stable  in  air,  but  the  solutions  are 
not  [11],     In  alkaline  solutions  [3,12]  the  salt  Is  converted 
into  pentacyanonltrltof errat e  (ll),  and  in  acid  solution 
hydrocyanic  acid  is  liberated  at  high  concentrations.  The 
solubility  of  the  salt  is  40  mg/lOO  ml  H^O  at  l6°C  [13]  and 
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57  mg/100  ml  at  43° C  [8],  The  salt  transforms  slowly  In 
direct  sunlight   [l4]  to  Fe^(CN)^g,  NO  and  Prussian  blue. 

The  specification  for  reagent  grade  sodium  nitroprusside 
is  as  follows 

Sodium  nitroprusside,  Na2[Fe( CN)^NO] •2H2O,  Formula 

wt.  297.97 

Insoluble  matter:     Not  more  than  O.i 
Chloride  (Cl"):     Not  more  than  0 
Sulfate  (SO^):     To  pass  test  (limit  about  0.01^). 
Microchemical  analysis  on  three  single  crystals  taken  at 
random  was  performed  at  NBS;,  with  the  results  -[15]: 

Iron  18.74  +  0.05^^ 

Carbon  20.2    t  0 ,3% 

Hydrogen        1.4    t  0 .1% 

Nitrogen      28.0    +  0.3^. 
The  analyses  confirm  that  the  theoretical  composition 
may  be  used  in  any  calculations  required  with  this  material. 

C .     Physical  Description 

The  salt  crystallizes  in  deep  red  orthorhombic  crystals, 

12 

belonging  to  the  space  group  '^^yT  "^^^  cell  constants, 

for  four  formula  weights  per  unit  cell,  are  [15]: 

a  =    6.17  +  0.03  A 

b  =  11.84  +  0.06  A 

c  =  15.43  -  0.08  A. 
The  x-ray  density  [I7]  is  I.76  g/cm3,   a  figure  which  is  in 
agreement  with  the  macroscopic  determination  [l8j  of  1.72  't 
0.02  g/cm^. 

The  faces  are  the  (110)  primary  rhombic  prism,  the  (010) 
primary  brachy-pinakoid  and  the  (Oil)  primary  brachydomal 
prism.     The  platelets  are  cut  along  the  100  plane  as  shown  in 
figure  2  and  3j   in  the  x-ray  st ereographic  projection  shown 
in  figure  4. 


■5f 

Unless  otherwise  specified,  all  error  limits  are  standard 
deviations  of  a  single  observation. 
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D.     Electronic  Structure 

The  bonding  In  the  nltroprusslde  anion  has  been  described 
by  Manoharan  and  Gray  [19]  In  terms  of  molecular  orbital  theory 
with  the  configuration  3d^*^^^  4s  ^'^"^^^p  0.251^     rj.^^  44.9^  4s 
character  differs  with  the  60fo  4s  character  obtained  from  the 
Walker,  Werthelm,  Jaccarlno  plot   [20].     This  discrepancy  Is 
under  current  Investigation  In  several  laboratories  [20,21,22], 

The  quadrupole  splitting  Is  consistent  with  the  assignment 
of  a  low  spin  configuration  for  Fe(ll).     From  the  asymmetry 

of  the  Intensity  of  the  2  transitions,  ±3/2      ±1/2,  ±1/2  ±1/2, 
the  electric  field  gradient  Is  axlally  symmetric  and  positive 
In  sign  [23]   (see  Appendix  B) . 


Figure  4.     X-ray  stereographlc  projection  of  sodium  nltro- 
prusslde single  crystal  along  the  001  axis. 
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3.     MOSSBAUER  SPECTROSCOPY  MEASUREMENTS 

A.     Description  of  Equipment 

For  the  calibration  of  the  standard  two  types  of  Mossbauer 
spectrometers  were  used:     the  tandem  spectrometer ^  using  a 
multiplex  technique  for  Int er-comparlson  of  the  platelets,  and 
a  high  accuracy  spectrometer  using  an  optical  Interferometer 
for  the  velocity  measurements. 

The  tandem  spectrometer  operates  In  the  constant  accelera- 
tion mode.     The  Doppler  velocity  Is  generated  by  an  electro- 
mechanical transducer,  using  a  loudspeaker  to  produce  the  motion, 
that  Is  mechanically  coupled  to  a  velocity  sensor  Incorporated 
In  a  feedback  loop.     The  spectrum  is  accumulated  In  a  multi- 
channel analyzer  operating  in  a  multlscaler  mode  [24].  The 
analog  voltage  signal,  varying  linearly  with  time  in  a  saw 
tooth  fashion.  Is  taken  from  the  address  scaler  to  provide  the 
input  signal  for  the  transducer.     This  method  guarantees  syn- 
chronization of  the  Doppler  drive  velocity  with  the  channel 
number  of  the  analyzer.     In  order  to  improve  the  accuracy  of 
the  spectrometer,  two  spectra  are  accumulated  simultaneously 
(multiplex).     The  Doppler  velocity  that  is  used  for  both  source- 
absorber  combinations  is  provided  by  one  electromechanical  drive 
unit   [25].     For  the  calibration,  the  standard  platelets  were 
compared  with  the  primary  standard  in  the  multiplex  system. 
Cobalt-57  in  a  two  mil  palladium  matrix  (0.25  in  diameter)  was 
used  for  the  source.     The  l4.4  keV  radiation  was  detected  by  a 
gas  flow  proportional  counter,  using  a  10^  methane,  SOfo  argon 
(P-IO)  gas  mixture^     A  moving  source  geometry  was  used,  with  a 
0.125  in  aperature,  in  order  to  reduce  any  error  due  to  the 
cosine  effect  [2?]. 

The  systematic  error  with  respect  to  velocity  of  the  above- 
described  spectrometer  is  usually  not  within  the  attainable 
limits  of  random  error.   Since  it  is  desirable  to  use  the  stand- 
ard for  calibration  of  Doppler  velocity  as  well  as  for  chemical 
shift,  a  high-accuracy  spectrometer  was  designed  to  measure 
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the  quadrupole  splitting  of  sodium  nitroprusside .     An  inter- 
ferometer can  be  used  with  minor  modifications  to  measure  the 
velocity  directly  in  terms  of  the  primary  standards  of  length 
(green  line  of  "^^^Hg)  and  time  (WW  NBS  standard  frequency), 
since  velocity  is  determined  by  a  measurement  of  distance  and 
time . 

The  interferometer  consists  of  a  stationary  light  source, 
beam  splitter,  reflector,   a  movable  second  reflector,   and  a 
photomultiplier  (see  figure  5).     The  light  Intensity  produced 
by  the  interferometer  varies  sinusoidally  with  the  optical 
path  difference  (AS)  between  the  two  reflectors-. 


Figure  5.     Optical  system  for  accurate  measurement 
of  velocity. 
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The  current  produced  by  the  photomultiplier  Is  represent- 
ed by  the  equation: 

I(AS)  =        +  Iq  cos  ^^^^  ,  where  (1) 

=  dark  current 

Iq  =  signal  current. 

If  the  second  reflector  is  not  stationary^  but  moving  with  a 
velocity  v,  the  photomultiplier  current  will  be  time  dependent: 
I(t)  =  Ip  +  Iq  cos  27r  (-^H  ,  where  (2) 
.  2v 

frequency  =         .  ^ 

For  a  velocity  of  1  um/s,  using  the        Hg  green  line  (X.  = 
5460.7^  a),  this  frequency  corresponds  to  3663.5  hertz.  This 
interferometer  thus  allows  velocity  to  be  determined  by  a 
frequency  measurement.     To  provide  the  feed-back  signal  for 
the  electromechanical  transducer,  the  fringe  count  frequency 
is  compared  with  the  driving  frequency  of  a  General  Radio 
Company  model  1162A.  synthesizer*  (standardized  against  WW)  by 
means  of  a  one  to  one  frequency  correlator  [28].     With  this 
arrangement,  the  velocity  can  be  controlled  and  measured  with 
high  precision  and  accuracy.     In  this  constant  velocity  mode, 
the  spectrum  of  a  primary  standard  was  obtained  by  varying  the 
frequency  synthesizer  through  the  range  of  0  to  5  khertz  in 
100  hertz  steps.     Gating  circuits  were  used  to  control  the 
fringe  and  standard  frequency  counters  for  positive  and  nega- 
tive velocities  and  to  stop  the  radioactivity  counting  scalers 
during  the  spectrometer  reversing  interval.     The  detector  sys- 
tem is  identical  to  the  one  used  in  the  constant  acceleration 
mode.     More  detailed  descriptions  of  this  spectrometer  appear 
in  reference  [29], 


The  identification  of  commercial  equipment  in  this  paper  does 
not  imply  recommendation  or  endorsement  by  the  National  Bureau 
of  Standards;  nor  does  it  imply  that  the  equipment  identified 
is  necessarily  the  best  available  for  the  purpose. 
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B.     Measurement  of  Mossbauer  Spectrum  Parameters 

From  a  lot  of  500  platelets,  60  were  randomly  selected, 
and  of  the  60,  ten  were  chosen  as  primary  standards.     The  ten 
platelets  were  determined  to  be  the  best  of  the  60  by  visual 
examination.     After  calibration  of  the  primary  standards  with 
the  optical  spectrometer,  one  of  the  ten  was  selected  for  use 
as  a  comparator  in  the  calibration  of  the  50  sample  platelets. 
The  platelet  selected  was  determined  to  be  the  one  closest  to 
the  average  of  the  quadrupole  splitting  of  the  ten  primary 
standards . 

The  spectrum  parameters  of  chemical  shift,  quadrupole 
splitting,  percent  effect  and  linewidth  were  measured.     All  of 
these  parameters  were  calculated  with  the  mathematical  tech- 
nique described  in  Appeddix  A. 

1.  Chemical  Shift 

The  chemical  shift  of  the  sodium  nitroprusside  platelet 
measured  with  the  Incident  radiation  perpendicular  to  the  100 
crystal  plane,   at  25.0°C  is  zero  cm/s  by  definition,   and  is 
determined  by  the  center  of  the  doublet  described  below.  The 
shift  of  an  iron  compound  measured  relative  to  sodium  nitro- 
prusside is  called  the  differential  chemical  shift  5^,  while 

designates  the  chemical  shift  for  the  standard.  The 

SL  Q 

relationship  for  calculating  differential  chemical  shift  is 

6     =6     -  6  ,  , 
o        X  std 

where  6     is  the  chemical  shift  of  the  unknown  compound.  The 
temperature  dependence  of  the  chemical  shift  is  3  x  10  cm 
s""^  °K"^  [3O]. 

2.  Electric  Quadrupole  Splitting 

The  nitrous  ligand  in  sodium  nitroprusside  produces  an 
electric  field  gradient  at  the  iron  nucleus,  which  results  in 
the  observed  doublet  spectrum  [31].     The  average  quadrupole 
splitting,  AE„,  for  the  ten  primary  standard  platelets  has  been 
determined  to  be  O.I726  ±  0.0002  cm/s  at  25.0° C  with  a  tempera- 
ture dependence  [30]  of  10"^  cm  s  °K 

11 


For  the  100  crystal  plane^  the  doublet  Is  symmetric 
energy  [32]   (Appendix  B) .     A  misalignment  of  2  degrees  in 
the  100  crystal  plane  of  the  standard  from  a  perpendicular  to 
the  gamma  ray  produces  an  asymmetry  of  S  =0.1. 

3.  Llnewidth 

The  linewidth  (r)  measured  for  the  standard  is  0,0305 

-  0.0004  cm/s.     The  effective  absorber  thickness  for  the 

14.4  keV  gamma  radiation  [33,34]  is  T^  =  3.77  mg/cm^. 

a 

Corrections  for  cosine  effect  [32]  were  made.  Using 
Margulies'  method  [33]  the  thick  absorber  line  broadening  is 
1.51,  yielding  a  corrected  line  width  of  0.0202  -  0.0004  cm/s 
(Appendix  E).     For  powered  sodium  nitroprusside  with  absorber 
thickness  of  25  mg/cm    iron,  the  halfwidth  is  5^  larger  than 
the  theoretical  width. 

4.  Magnitude  of  Resonance  Effect 

The  standard  absorber  exhibits  a  percent  effect  (£)  of 

(TV 

17  ±  2fo  using  10  mCi  -^'Co  on  a  Pd  matrix,  with  an  argon-methane 
(P-IO)  gas  flow  proportional  counter.     The  single  channel 
analyzer  was  set  to  pass  9Q%  of  the  l4.4  keV  radiation. 

5 .  Gamma  Ray  Attenuation 

The  attenuation  of  the  crystal  for  l4.4  keV  (0.864  A) 
radiation  is  Q^fo,  as  determined  using  standard  y-spectrometric 
techniques.     Prom  Appendix  E  the  attenuation  can  be  seen  to  be 
due  to  the  iron  mass  absorption. 

C.  Computations 

Mossbauer  spectra  are  particularly  suited  to  resolutipn 
by  mathematical  methods.     Machine  digital  computation  is  being 
used,  hence  errors  in  data  manipulation  are  minimized  and 
comprehensive  random  error  estimates  can  be  made  of  the  spec- 
trum parameters.     The  mathematical  techniques  used  for  resolv- 
ing the  Mossbauer  spectra  are  described  in  Appendix  A.  An 
example  of  the  output  of  the  computer  program  called  PARLORS 
IV  is  given  in  Appendix  C. 
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4.     ANALYSIS  OF  RANDOM  AND  SYSTEMATIC  ERRORS 

Since  the  chemical  shift  of  this  standard  reference 
material  Is  defined  to  be  zero  there  Is  no  systematic  error 
associated  with  It,    Therefore,  systematic  errors  are  associ- 
ated only  with  the  quadrupole  splitting.     Since  we  define  the 
midpoint  between  two  peaks  of  the  spectrum  of  the  standard  as 
zero  differential  chemical  shift ,  the  random  error  in  the  quad- 
rupole splitting  reflects  that  associated  with  the  measurement 
of  the  zero  point  for  differential  chemical  shift.     Ten  primary 
standards  were  selected  to  be  calibrated  on  the  high  accuracy 
optical  spectrometer.    The  average  quadrupole  splitting  of 
these  ten  results  Is  the  value  reported  on  the  certificate. 
The  secondary  standards  were  then  calibrated  from  one  of  the 
primary  standards  with  the  use  of  the  tandem  spectrometer  to 
determine  the  random  error  of  both  the  quadrupole  splitting 
and  the  zero  point  of  differential  chemical  shift  for  the  50 , 
samples  that  represent  the  lot  of  500  platelets.    The  random 
errors  In  the  chemical  shift  and  quadrupole  splitting  are  the 
standard  deviations  of  a  single  determination. 

A.     Random  Errors  In  the  Calibration  Procedure 

Random  errors  resulted  from  the  comparison  technique  and 
the  optical  spectrometer  measurements.    The  random  errors  In 
the  latter  will  be  analyzed  first. 

For  velocity  measurements  using  the  high  accuracy  spectro- 
meter, the  distance  that  the  absorber  traversed  in  a  unit  of 

time  was  determined.     Since  the  Mossbauer  transition  takes 

-7 

place  In  a  time  Interval  of  the  order  of  10     s,  the  unit  of 
time  for  this  measurement  should  be  of  this  magnitude.  An 
Interf erometrlc  technique  using  visible  light  reduces  the 
measurement  of  velocity  to  a  frequency  measurement.     For  this 
method  the  unit  of  length  is  the  "^^^Hg  green  line  (y  =  5460.74 
±  0.01  A  corrected  for  air)   [32].     By  comparing  the  NBS  stan- 
dard frequency  (accuracy  of  1  part  in  10^)  with  the  fringe- 
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count  frequency  generated  by  the  interferometer  using  the 
198 

Hg  green  line  the  velocity  can  be  measured.  To  evaluate 
the  accuracy^  the  relationship  between  velocity  and  fringe- 
count  frequency  (f)  is 

v=^  (1) 

since  the  optical  path  is  twice  the  traverse  of  the  spectro- 
meter.    The  fractional  error  in  velocity,  E(v),  is  given  by: 

E(v)  =  [E{\f  +  E{ff]K  vJhere  (2) 

E(\)  =  fractional  error  in  wave 
length  =  2  X  10-6 

E(f)  =  fractional  error  in  frequency. 

The  interferometer  is  capable  of  measurement  to  within 
1/2  fringe,  therefore  the  fractional  error  in  the  frequency 
is  given  by: 

E(f)  =^  =  i  (^)  =^  (3) 

substituting  into  equation  2 

E(v)  =   [4  X  10-12  ^  (^)2ji  .^^ 

for 

V  =  1  cm/s 

E(v)  =  2  X  10"^ 

Equation  (4)  is  an  expression  for  the  limiting  random  error. 
This  expression  also  indicates  that  the  interf erometric  method 
will  not  work  at  zero  velocity  because  the  error  becomes 
infinite. 

The  above  consideration  does  not  take  into  account  the 
limitations  on  the  traverse  of  the  spectrometer,  due  to  visible 
light  coherent  path  length  and  the  Inverse  square  law  of  radia- 
tion losses  if  long  traverses  are  used. 

Therefore,  a  more  realistic  estimate  of  the  error  can  be 
obtained  from  the  error  in  the  fringe  counts.     Since  the 
electronic  circuitry  is  capable  only  of  counting  to  the  nearest 
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fringe,  the  error  In  velocity  Is  then  dependent  upon  the  num- 
ber of  fringe  pulses  (f)  per  centimeter  of  travel  (d). 

^  =  (\/2Hl/2}  (5) 

The  factor  l/2  arises  from  the  fact  that  the  electronics 
count  half-cycles.     For  the  one  centimeter  distance  of  travel 
used  In  this  work,  equation  (5)  yields  8  x  10    fringe  pulses. 
Hence  the  random  error  In  velocity  for  one  centimeter  travel 
Is  one  fringe  pulse. 

The  above  considerations  assume  the  velocity  correction 
circuitry  to  be  capable  of  Instant  correction  when  a  velocity 
error  Is  sensed.     However,  the  mass  of  the  moving  system  re- 
duces the  correction  time,  which  produces  a  random  error  If 
the  velocity  corrections  are  random,  and  a  systematic  error  If 
the  correction  signs  are  largely  positive  or  negative.  This 
systematic  error  will  be  discussed  below.     The  random  error 
could  only  Increase  the  line  width  determination,  and  the 
amount  of  Increase  In  line  width  does  not  significantly  in- 
crease the  random  error  in  the  quadrupole  splitting  determina- 
tion 

Since  the  tandem  system  was  used,  only  the  differential 
linearity  (e.g., the  line  width  divided  by  the  separation  of 
the  doublet)  is  important.     This  differential  linearity  was 
found  to  contribute  an  insignificant  amount  to  the  total  ran- 
dom error  for  the  calibration  measurements  [24], 

1.     Optical  Spectrometer  Drift 

The  optical  spectrometer  drift  would  contribute  the 
largest  error  in  the  calibration  measurements,   since  this 
would  be  directly  reflected  in  the  random  error  of  the  chemical 
shift.     Evaluation  of  the  magnitude  of  this  error  by  measuring 
one  of  the  primary  standards  over  a  one  month  period  indicated 
a  standard  deviation  of  the  chemical  shift  of  ±  0.02  cm/s  from 
50  measurements.     However,  by  using  the  tandem  spectrometer 
this  error  due  to  drift  was  reduced  to  a  negligible  amount. 

15 


2 ,     Counting  Statistics 

Prellmlnapy  measurements  on  the  sodium  nltroprusslde 
crystals,  and  the  economics  of  data  accumulation  time  dictated 
that  a  maximum  of  ^  x  10^  counts  per  channel  (N^)  for  200 
channels  be  accumulated  for  each  platelet.     This  would  corre- 
spond to  a  standard  deviation  of  a  single  determination  of 
10      cm/s  due  to  counting  statistics.    This  was  further  con- 
firmed by  the  residuals  obtained  from  a  least  squares  fitting 
of  the  experimental  data  points  to  a  Lorentzlan  line  profile, 
discussed  In  Appendix  A. 

B.     Systematic  Error  In  the  Determination  of  the  Electric 
Quadrupole  Splitting 

The  error  In  the  velocity  measurement  Is  limited  by  the 
Instrumentation,  due  to  the  resolution  of  the  one  to  one 
frequency  comparator  and  the  time  response  of  the  system.  The 
open  loop  gain  as  a  function  of  frequency  decreases  6  dB  per 
octave  and  the  slgnal-to-nolse  ratio  limits  the  response  to 
300  hertz.     The  frequency  correlator  response  time  Is  Inversely 
proportional  to  the  frequency,  but  the  response  time  can  be 
Increased  by  an  Integrating  network  at  the  output  of  the  com- 
parator so  as  to  match  the  decrease  In  response  of  the  drive 
unit . 

The  error  due  to  the  system  response  could  be  measured 
directly  from  the  comparison  (beat)  of  the  fringe  frequency 
with  the  reference  frequency,  using  an  oscilloscope  or  by 
counting  the  y  radiation  at  a  velocity  of  0.5  T  away  from 
resonance  as  a  function  of  the  velocity  correction  signal. 
These  random  fluctuations  never  exceeded  10^  of  the  reference 
frequency,  which  would  correspond  to  a  maximum  error  of  0.1  x 
1/4X  •  f  or  typically  2fo  of  the  line  width  of  the  standard  In 
the  worst  case.     Since  the  best  frequency  was  random,  this 
error  could  only  Increase  the  line  width  determination,  and 
not  the  separation  of  the  peaks.     The  systematic  error  as  de- 
termined from  the  total  fringe  count  with  respect  to  the 
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reference  frequency  counts  was  found  to  be  within  the  standard 
deviation  of  ±  5  in  10^  fringes. 

The  cosine  effect  was  minimized  by  using  the  proper 
aperture.     For  the  measurements  on  the  primary  standards,  a 
0,25"  diameter  source  and  detector,  separated  by  2.0"  were 
used,  and  the  cosine  effect  correction  was  made  for  this 
geometry.    The  alignment  of  source  and  absorber  gave  negligible 
error , 

In  order  to  determine  any  long  term  stability  effects  on 
the  standard,  a  series  of  spectra  were  taken  with  the  high 
accuracy  optical  Mossbauer  Spectrometer  over  a  period  of  one 
year.    No  detectable  changes  in  any  of  the  parameters  were 
observedo     However,  it  is  always  considered  advisable  to  store 
the  crystal  in  a  dark  place  away  from  any  sources  of  extreme 
heat,  cold  or  humidity. 

C,  Summary 

A  systematic  study  of  the  sources  of  random  errors  indi- 
cates that  this  error  is  represented  entirely  by  the  limitation 
in  the  total  number  of  Mossbauer  gamma  rays  detected.     A  lower 
random  error  can  be  obtained  at  the  expense  of  increased  mea- 
suring time.    The  limiting  random  error  would  be  that  caused 
by  spectrometer  drift.     A  similar  study  of  systematic  errors 
indicates  that  the  accuracy  is  well  within  the  limits  imposed 
by  the  random  error;     consequently  it  can  be  said  that  if  any 
systematic  error  exists  it  is  less  than  0.0002  cm/s  for  the 
quadrupole  splitting., 
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5.     THE  USE  OF  THE  STAITOARD 

Procedure  for  Measuring  Chemical  Shift  and  Reporting 
Differential  Chemical  Shift 

The  standard  crystal  can  be  used  to  calibrate  the  spec- 
trometer before  or  after  a  measurement,  or  the  standard  spec- 
trum and  the  spectrum  to  be  calibrated  can  be  accumulated 
simultaneously  If  a  tandem  spectrometer  Is  available.  Since 
the  unknown  Is  measured  relative  to  the  standard,  the  source 
parameters  do  not  affect  the  determination  of  the  differential 
chemical  shift.     For  comparison,  the  various  sources  and  their 
respective  chemical  shifts  are  shown  In  figure  1. 

For  detailed  directions  on  the  use  of  this  standard,  see 
the  Certificate  of  Calibration  on  page  22. 
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6 .  SUMMARY 

The  characteristics  of  sodium  nitroprusside  fully  meet 
the  requirements  for  an  acceptable  Mossbauer  standard.  Its 
chemical  stability,  uniformity  and  doublet  spectrum  give  the 
required  accuracy  needed  for    this  spectroscopy.     The  re- 
sults of  this  study  are  best  summarized  in  the  Certificate  of 
Calibration  which  follows  this  page. 

The  use  of  this  standard  and  the  use  of  consistent  re- 
porting procedures  of  the  spectrum  parameters  should  provide 
stability  to  the  growing  field  of  Mossbauer  spectroscopy.  In 
addition,  it  should  accelerate  efforts  to  quantitatively  unde: 
stand  the  spectrum  parameters  by  providing  the  theoretician 
with  highly  reliable  data. 

This  Standard  Reference  Material  can  be  purchased  from 

Office  of  Standard  Reference  Materials 
National  Bureau  of  Standards 
Washington,  D.  C.  20234 

Attention: 

Standard  Reference  Materials  Unit. 
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Directo,   (^j^j^iifi^^i^  QIaltbrattnn 

Standard  Reference  Material  725 
for 

Mossbauer  Differential  Chemical  Shift  for  Iron-57 
Sodium  Pentacyanonitrosylferrate(II)  Dihydrate* 

(Sodium  Nitroprusside) 

This  Standard  Reference  Material  was  prepared  from  a  single  crystal  of  the  compound  sodium  pentacyano- 
nitro9ylferrate( II)  dihydrate ,  Na[Fe(CN)5N0]  •  2H2O.  The  purity  of  this  compound,  also  known  as  sodium 
nitroprusside,  meets  the  specifications  of  the  American  Chemical  Society  for  reagent-grade  materials  as 
verified  by  quantitative  analysis  of  the  main  constituents  ,  but  should  not  be  considered  as  entirely  free  from 
impurities  such  as  heavy  metals.  It  is  in  the  form  of  a  platelet  of  dimensions  1  x  1  x  0.0775  ±  0.003**  cm 
that  has  been  cut  from  a  large  single  crystal.  The  1  x  1  cm  surface  is  parallel  to  the  100  crystal  plane  within 
±  2  degrees  of  arc.  The  opposite  1  x  1  cm  surfaces  are  parallel  to  within  0.001  cm  with  a  surface  finish  of  20 
microns.  The  natural  iron  concentration  is  25.0  ±  1  mg/cm2.  The  chemical  shift  of  this  Standard  is  compared 
to  that  of  the  National  Bureau  of  Standards  Primary  Standard  which  has  the  value  of  zero  cm/s  differential 
chemical  shift.  This  Standard  has  an  average  value  for  the  chemical  shift  of  0.0000  ±  0.0002  cm/s  at  25.0°C. 
This  uncertainty  is  expressed  as  the  standard  deviation  of  a  single  determination  derived  from  single  measure- 
ments on  50  platelets.  The  National  Bureau  of  Standards  dual-spectra  Mossbauer  spectrometer  [l,2]  was 
used.  The  resonant  spectra  have  considerable  line-broadening  due  to  the  thickness  of  the  absorber.  The 
experimental  line  width  (full  width  at  half  maximum)  is  0.0305  ±  0.0004  cm/s.  The  line  width  corrected  for 
thickness  broadening  is  0.0202  cm/s.  The  Mossbauer  effect  is  17  percent. 

To  allow  the  user  to  calibrate  the  velocity  scale  of  a  Mossbauer  spectrometer,  the  electric  quadrupole 
splitting  of  the  Standard  Reference  Material  was  measured.  The  velocity  scale  of  the  dua'-spectra  comparison 
spectrometer  was  calibrated  by  measurement  of  the  electric  quadrupole  splitting  of  the  Primary  Standards 
using  an  optical  interferometric  technique.  The  average  value  of  the  electric  quadrupole  splitting  for  this 
Standard  is  0.1726  ±  0.0002  cm/s  at  25.0°C  where  the  uncertainty  is  expressed  as  the  standard  deviation  of  a 
single  determination  derived  from  single  measurements  on  10  platelets. 

The  details  of  the  instrumentation  and  of  preparation  and  measurement  of  the  sodium  nitroprusside 
crystals  are  given  in  NBS  Miscellaneous  Publication  260-13  entitled  "Mossbauer  Spectroscopy  Standard  for 
Chemical  Shift  of  Iron  Compounds." 

The  single  crystal  platelets  of  sodium  nitroprusside  were  prepared  by  the  Isomet  Corporation  of  Palisades 
Park,  New  Jersey.  The  calibrations  of  the  crystals  were  made  within  the  National  Bureau  of  Standards 
Institute  for  Materials  Research,  by  J.  J.  Spijkerman,  F,  C.  Ruegg,  D.  K.  Snediker,  and  W.  L.  O'Neal  of  the 
Radiochemical  Analysis  Section ,  James  R.  DeVoe,  Chief. 

*  This  nomenclature  is  based  upon  the  structure  of  this  compound  as  determined  by  Mossbauer  spectroscopy. 

**  Unless  otherwise  specified  all  uncertainties  are  standard  deviation  of  a  single  determination. 

References: 

[1]  Rad  iochemical  Analysis:  Activation  Analysis,  Instrumentation,  Radiation  Techniques ,  and  Radioisotope 
Techniques,  July  1963  to  June  1964,  James  R.  DeVoe,  Editor,  NBS  Tech.  Note  248 ,  pp  25-37  (1964).  Avail- 
able from  Superintendent  of  Documents ,  Government  Printing  Office,  Washington,  D.C.  20402,  50  cents  a  copy. 

[2]  Radiochemical  Analysis:  Activation  Analysis,  Instrumentation,  Radiation  Techniques ,  and  Radioisotope 
Techniques,  July  1964  to  June  1965 ,  James  R.  DeVoe,  Editor ,  NBS  Tech.  Note  276  ,  pp  74-1 10  ( 1966).  Avail- 
able from  Superintendent  of  Documents,  Goviemment  Printing  Office,  Washington ,  D.C.  20402,  $1.00  a  copy. 


Washington,  d.  c.  W.  Wayne  Meinke,  Chief 

January  30,  1967  Office  of  Standard  Reference  Materials 

(This  certificate  supersedes  certificate  of  4-6-66) 
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Directions  for  Use 


The  Standard  Reference  Material  in  the  form  of  a  single  crystal  platelet  is  sandwiched  between  two 
pieces  of  a  4  mil  polyethylene  film.  Due  to  the  advisability  of  keeping  it  free  from  moisture  it  is  recom- 
mended that  the  platelet  not  be  removed  from  this  container.  A  suitable  mounting  consists  of  two  concen- 
tric aluminum  rings  approximately  1.25  inches  in  diameter.  The  crystal,  encased  in  the  plastic  envelope, 
is  placed  between  the  rings  and  the  assembly  is  fastened  together  with  small  screws,  therelDy  clamping  the 
envelope  firmly  between  the  rings.  Both  rings  have  an  inside  diameter  of  0.75  inch  to  expose  the  crystal  and 
to  prevent  crushing  of  its  edges.  This  was  the  procedure  used  in  taking  all  of  the  data  at  NBS  on  this 
Standard  Reference  Material. 

Place  the  Standard  Reference  Material  as  the  absorber  at  25.0°  C  in  the  Mossbauer  spectrometer  and 
take  sufficient  transmitted  counts  of  the  14.4  kev  gamma-ray  from  an  iron-57  source  to  obtain  counting 
statistics  to  the  required  degree  of  uncertainty.  Determine  the  peak  position  of  the  electric  quadrupole 
splitting.  The  use  of  digital  computation  techniques  are  recommended.  Divide  the  distance  between  the 
peaks  by  two  and  assign  this  as  the  value  for  chemical  shift  indicated  above.  Replace  the  Standard  crystal 
with  the  compound  whose  chemical  shift  is  to  be  measured  as  the  absorber  in  the  spectrometer.  Measure 
the  peak  position  parameters  as  a  difference  between  that  of  the  unknown  compound  and  the  mid-point  be- 
tween the  peaks  of  the  Standard. 
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APPENDIX  A 

COMPUTER  CURVE  FITTING  PGR  MOSSBAUER  SPECTRA 


The  mathematical  analysis  described,  below  resolves  a 
series  of  Mossbauer  peaks  superimposed  on  a  baseline  parabola. 
Lorentlzlan  functions  are  used  for  the  Mossbauer  absorption 
profile,  and  since  the  equations  Involved  In  finding  the  best 
fit  are  non-linear,  a  solution  cannot  be  obtained.  However, 
by  making  several  approximations  and  substitutions,  linear 
equations  result  which  can  be  solved  for  approximate  results. 
The  nature  of  these  approximations  requires  initial  estimates 
of  peak  position  and  half  width  (half  width  at  half  maximum) 
to  be  made  manually.     Once  the  equations  are  solved,  results 
of  the  first  solution  can  be  used  as  new  input  estimates,  and 
the  solution  procedure  repeated.     By  iterating  in  this  fashion, 
a  solution  is  obtained. 

The  basic  form  of  a  curve  consisting  of  n  Lorentzlans  and 
a  parabola  is 


where  the  first  bracketed  term  is  the  sum  of  the  n  Lorentzlans 

and  the  second  term  the  equation  for  the  parabola.  Here  A^,  P^^ 
and  B.   are  the  amplitude,  position,  and  half  width  of  the  ith 
(Lorentlzlan)^  and  E,  F  and  G  have  their  usual  meaning  .  These 
3n+3  parameters  are  the  unknowns  to  be  solved  for  by  the  pro- 
gram.    Since  a  successive-approximation  method  is  used,  input  to 
the  computer  must  Include  a  convergence  criterion,   specifying  the 
lower  limit  parameter  corrections  desired  in  the  solution^  and 


This  method  is  used  in  a  series  of  computer  programs  called 
PARLORS  written  in  FORTRAN  4.1  for  the  CDC  36OO  [25]. 
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a  maximum  number  of  iterations  (e.g.  a  limit  on  the  number  of 
times  the  program  will  iterate  in  the  event  of  diverging  solu- 
tions) as  well  as  the  required  estimates. 

Once  a  suitable  number  of  iterations  is  obtained  the 
program  prepares: 

1)  a  subtract  matrix,,  containing  distance  between  any  two 
peaks  above  the  diagonal,  peak  position  along  diagonal  and  the 
centroid  of  any  two  peaks  below  the  diagonal; 

2)  an  area  fraction  matrix,   containing  the  ratio  of  the 
area  under  two  peaks  off  the  diagonal,   and  the  fraction  of  the 
area  under  one  peak  to  the  total  area  on  the  diagonal; 

3)  an  error  analysis;  and 

4)  several  types  of  plots. 

METHOD  OF  SOLUTION 

To  find  the  best  least  square  fit  to  a  set  of  data  (x  , 
.  J 
y.)  the  quantity 
J 

k 

^  =    ^  [^j(obs)-yj(theor)]^ 
J-^1 

must  be  minimized.     (For  most  sets  of  data,  k  is  equal  to  200 
channels.)  From  (1)  the  expression  for  Y^(^-^YieoT)  substituted: 

k 


0=1 


y 


j(obs)' 


1  + 


+  E  +  Fx 


Gx 


2 


To  linearize  the  equation  the  following  definitions  are  made 


h. 

1 


where 


« 


(5) 
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1 


(6) 


where        and        are  the  estimated  values  of  h^  and  P^,  the 
peak  height  and  position  with  error  6^  and  y^,  respectively, 
Using  (4) 

(7) 


1  + 


to  linearize  the  equations^  (5)  and  (6)  are  substituted  Into 
(7)  .  The  result  Is  expanded  by  the  binomial  theorem,  and  all 
higher  terms  (non  linear)  of  6^,  y^  are  neglectedo 

A.  A. 


1  + 


A.. 

(8) 


1  + 


1  + 


Let 


becomes 


Ij  "  ^l(^j"-^l^^^  ^^^^  right  side  of  equation  (8) 


Substituting  (8)  Into  (3) 
k  r  n 


S  = 


j=i  1+1 


(9) 
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Letting  =  2A^H^y^  and        =  A^6^ 


k 

J-1 


i=l  ^Ij  ij 


+  (10) 


.  .  0   E  +  FXj  +  GXj 

The  unknowns  are  now  A^fC^,D^,  for  1  =  l,2,...n  and  E,F,G  used 
in  the  minimization  of  S  results  in  a  set  of  3n+3  simultaneous 
linear  equations,  which  are  the  usual  normal  equations  of  least 
squares.     The  following  manipulations  simplify  setting  up  these 
equations,   and  in  fact  follow  the  method  used  in  the  computer 
program: 

Let  y'  be  the  row  vector  (prime  denotes  transpose) 

y'  =  (yi  yg  •••yk^  ^^^^ 

denoting  the  k  observations  and  let  Z'  be  a  row  vector 


Z  =  (A-L.C-^.D^.A^.C^oo.A^^C^^D^^E^P^G) 

denoting  the  3n+3  unknowns  and  let 


T  = 


1 

1 

(Xi-P2)^ 

2  ^ 

2  ' 

1 

X2-P2 

(Xg-Pg)^ 

«2l' 

2 

2  ^ 

«22 

2 
Q22 

1-    V^i  (VPi)""    1    V^2  (V^s)^ 


'kl 
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^k2  q: 


k2 


^k2 


(12) 


l^X^^X^ 

(13) 


So  that  the  normal  equations  become 

T'TZ  =  T'y 

29 


(14) 


The  matrix  equation  (l4)  Is  then  solved  for  the  values  of  Z. 
Then 


A.  =  Z 


31-2 


,    .  ^31-1 
^1  2H,Z3,_.2 


^31  ^31 


(15) 


i  A. 


'31-2 


3n+l 


F  =  Z 


3n+2 


G  =  Z 


3n+3 


From  these  corrections,   y.   and  6.^3  better  estimate  of  P.  and 


h.    (and  thus  B. )  are  obtained. 


1-    -  -1 

(Note  that  because  of  the  ap- 


proximations made  the  solutions  are  not  exact).     If  more  ac- 
curacy Is  desired,  reiterate  by  replacing  the  old  estimates 
with  the  corrected  estimates  and  repeat  the  whole  procedure. 

ERROR  ANALYSIS 
Once  the  value  of  each  unknown  has  been  obtained,  the 
theoretical  value  of  the  data  points  (y^.^heor^  calcu- 
lated.    The  mean  square  residual. 


St 

M-1  ^ 


^J-obs"  ^j-theor 


/(k-3n-3) 


(16) 


describes  the  precision  of  fit,     (k  Is  the  number  of  data 
points;  and  3n+3  the  number  of  parameters  solved  for.)  The 
varlance-covarlance  matrix  for  Z  Is  then 


Var(Z)  =  (T'T) 


■1  2 

a 


(17) 


The  variances  of  A^,  E,F,  and  G  are  given  by  the  diagonal 
terms  of  (17)  and  the  covarlance  of  Z.   and  Z.  by  the  l,j  terms 
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Using  propagation  of  error  formulae  and  (15) 


Var 


!!iy+    _  gc°v(Ai,Dj 


(18) 


Then  the  approximation  Var  (h. 


(6^)  Is  made  and  since  from 


=  1/  h^. 


Var(B, ) 


rdB."]^  Var  (6 


(19) 


Since  Yj_  =  Var(Yj_)  depends  on  Var  (H^).  Because 


h^  =  H^^  the  equation 


(20) 


Is  used  In  the  error  analysis  to  take  account  of  the  eri-or  In 
half  width. 

Because  A^^^  C^,  and  J)^  are  determined  from  the  same  data 
they  will  Involve  the  same  random  errors  and  hence  be  correlated. 
This  correlation  Is  expressed  by  the  terms  of  the  varlance- 
covarlance  matrix  as  follows: 


M  = 


■Var(Aj_)  Cov(A^C^) 
Cov(A^C^)  Var(C^) 


Cov(A^D^)  Cov(C^Dj_) 


The  variances  and  covarlances  of 


Cov(A^D^)  - 

cov(c^D^ ; 

Var(D^) 


(21) 


[0 


0] 


(22) 
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and 


2A^H^-2D^  =  [2Hj_       0  -2] 


are  given  by 


0 

1 

0 

0 

M 

1 

2H, 

0 

-2 

0 

-2 

(23) 


(24) 


Var(C^)     [2H^Cov(A.  ,C^)-2Cov(C^,D^)] 
[2H^Co V  (     ,     )  -2 Co V  (     , D .  )][4h .  ^Var  ( A^  )  -8h^ Gov  (  A^  , D  .  )  +  4Var ( D .  )] 


Substituting  the  appropriate  matrix  elements  into  equation 
(l8)  gives 

2, 


Var(Yj_)  =  (y^)' 


Var  (     )      H^^'Var  ( A^  )  -2H^ Gov  ( A^D^  )+Var  ( 
— +  o  


A.H.-D. 


(25) 


2[H^Cov(A^,C. )-Gov(D^,G. )] 
^l^VV^V 


Finally  the  approximation  Var  (Pj_)  =  Var  ( y  j_ )  Is  made.  By 
the  above  formulas  the  values  and  variances  of  all  the  para- 
meters in  (1)  can  be  determined. 
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APPENDIX  B 


ANGULAR  DISTRIBUTION  OF  RADIATION  IN  THE 
MOSSBAUER  EFFECT 


We  consider  the  particular  case  of  the  Mossbauer  transi- 
tions In  Pe-57.     These  transitions  occur  between  the  ground 
state  of  spin- 1/2  and  the  first  excited  state  (l4.4  kev)  of 
spin- 3/2. 

The  spln-3/2  state  possesses  4  magnetic  substates  (^3/2, 
^1/2)  so  that  there  are  6  possible  Mossbauer  transitions: 

9/4(l+cos^e) 


3  sln^e 


3/2 

-» 

1/2 

Am  =  -1 

-3/2 

-* 

-1/2 

=  +1 

1/2 

-¥ 

1/2 

=  0 

-1/2 

-* 

-1/2 

=  0 

-1/2 

-¥ 

1/2 

=  +1 

1/2 

-¥ 

-1/2 

=  -1 

3/4(l+cos^e) 


Here  Am  represents  the  z-component  of  angular  momentum  which 
Is  carried  off  by  the  photon,  or  gamma  ray.   In  the  transition. 

The  angular  distribution  of  radiation  Is  derived  with  the 
aid  of  the  theory  of  multiple  fields.     In  this  particular  ex- 
ample we  are  concerned  with  transitions  which  Involve  magnetic 
dlpole  radiation.     The  angular  distribution  Is  given  by  the 
equation  [36] 

Z,  (e,$)  =  X*    •  X, 

=  V2[l-f[^]|Y,^,,,|^l/2[l-f[^] 

I  Y  |2   ,  |2 
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In  this  equation  6  Is  the  angle  between  the  gamma  ray  and  the 
quantization  axis,   and  $  is  the  azimuthal  angle  about  this 
axis;   I  is  the  multipole  of  the  radiation  and  m  is  the  z-com- 
ponent  of  angular  momentum  of  this  radiation;  the  Y^^  (9^^) 
called  the  vector  spherical  harmonics  and  are  related  to  the 
spherical  harmonics  by  the  relation 

X   (e,$)  =  2  Y   (e,$)     S  =  iti  ?  X  V 

m  m 

where  L  is  the  quantum  mechanical  angular  momentum  operator. 
The  physical  significance  of  the  vector  spherical  harmonic 
X_^^(0,$)  is  that  it  completely  determines  the  angular  depen- 
dence of  the  electric  and  magnetic  multipole  fields. 

The  relative  transition  probabilities  of  the  various 
transitions  are  determined  by  the  Clebsch-Gordon  coefficients 
(also  called  vector  addition  coefficients,  Wigner  coefficients, 
or  C-coeff icients) .  The  probability  that  a  transition  will  occur 
is  given  by         , ( JM;m,m' )   (notation  of  Blatt-Weisskopf  [36]), 
In  this  notation  I  is  the  spin  of  the  photon  (=l)  and  is 
the  spin  of  the  ground  state;  J,  M  are  the  spin  and  z  component 
of  spin  for  the  excited  state;  m  is  the  z-component  of  angular 
momentum  carried  off  by  the  photon  in  the  transition  and  m'  is 
the  z-component  of  spin  for  the  ground  state.     Clearly,  M  =  m 
+  m'   (conservation  of  z-component  of  spin).     The  transitions 
and  their  relative  probabilities  in  the  cases  when  the  degen- 


eracy  of  the  ma^ 

^netic 

substat es 

is  resolved  and  unresolved 

indicated  below: 

Transition 

Relative  Probability 

1/2  3/2, 

Am  = 

+1, 

3 

-1/2  ->  1/2, 

Am  = 

+1, 

1 

1/2  1/2, 

Am  = 

0, 

2 

-1/2  -1/2, 

Am  = 

0. 

2 

1/2  -1/2, 

Am  = 

-1. 

1 

-1/2  -3/2, 

Am  = 

-1. 

3 
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Transition 


Relative  Ppobability 


+3/2  ^  ±1/2 
±1/2  4*  +1/2 


1 


1 


In  the  case  when  the  degeneracy  is  unresolved^  the  angular 
distribution  of  radiation  and  the  relative  transition  probabil- 
ities are  given  by: 


Now  in  the  case  of  single  crystals  the  absorption  spectrum 
will  possess  a  pronounced  asymmetry.     This  is  due  to  the  fact 
that  the  distribution  of  radiation  is  affected  by  the  crystal 
structure  in  a  complicated  way  ( Go I'danskii-Karyagin  Effect  [37], 
which  varies  for  different  crystals.     In  this  case  the  distri- 
butions given  above  must  be  multiplied  by  the  so-called  Lamb- 
Mossbauer  or  Debye-Waller  factor  f(0)  which  depends  on  the  angle 
6  as  well  as  on  the  temperature  and  the  energy  of  transition. 
It  is  clear  that  the  factor  f(0)  acts  as  a  kind  of  distribution 
function  which  selects  only  the  Mossbauer  transitions,  i.e., 
only  recoil-free  transitions. 

From  figure  6,  the  electric  field  gradient  vector  is  at  an 
angle  of  54°  with  the  100  crystal  plane.  The  -5/2  state  is  of 
higher  energy  that  the  -1/2  state  and  the  electric  field  gradi- 
ent is  positive.  For  the  001  and  010  planes,  the  intensity 
ratio  is  1.55  and  0.65,  respectively.  The  calculated  values  of 
figure  6  differ  from  the  experimentally  observed  values  of  1.55 
and  0.74,   indicating  the  possibility  of  a  Gol'danskii-Karyagin 


±3/2  ±1/2 
+1/2  +1/2 


3/2(l+cos^e) 
1+3/2  sin^e 


1 


1 


Effect. 
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Figure  60     Peak  intensity  assymmetry  for  the  -3/2->-l/2 
doublet  in  iron  compounds  as  a  function  of 
the  angle  between  the  electric  field  gradient 
and  the  incident  radiation. 
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APPENDIX  C 

FOPIMAT  FOR  PRESENTATION  OF  EXPERIMENTAL  MOSSBAUER 
SPECTRA  AND  PARAMETERS 

Typical  of  every  new  spectroscopic  discipline,  there  has 
been  a  variety  of  presentations  of  spectra  using  different^ 
nomenclature  and  units  for  the  spectrum  parameters.     For  pur- 
poses of  encouraging  spectroscopist s  to  report  their  data  in 
a  uniform  manner  a  suggested  format  and  nomenclature  are  de- 
scribed below.     It  is  hoped  that  these  reporting  procedures 
■will  serve  as  a  basis  for  those  accepted  by  pertinent  govern- 
ing bodies  such  as  the  International  Union  of  Pure  and  Applied 
Chemistry. 

Nomenclature  for  the  pertinent  Mossbauer  parameters  are 
presented  in  Table  1.     Certain  explanations  are  necessary. 
The  electric  quadrupole  splitting  can  best  be  represented  by 
subscripting  the  magnetic  quantum  numbers  of  the  transitions. 
Since  the  separation  between  two  peaks  must  be  represented  as 
the  difference  between  two  sets  of  transition  energies,  the 
value  of  AEq  must  be  followed  by  a  subscripted  matrix  indicat- 
ing the  two  transitions.     The  magnetic  quantum  numbers  in  the 
excited  and  ground  state  of  transition  1,   (rn^jm^)  and  that  of 
transition  2,(m^,m2)  uniquely  define  AE^. 

The  line  width  r  obtained  experimentally  is  to  be  con- 
trasted with  r    the  theoretical  line  width.     Corrections  for 
o 

thickness  of  the  absorber  [33]  for  the  cosine  effect  [27]  and 
other  effects  Including  instrument  broadening  may  be  applied 
to   r  in  order  to  approach  r^. 

Reporting  the  magnitude  of  the  effect  requires  great  care 
if  laboratory  Independency  on  measurement  is  to  be  retained. 
The  expression  given  requires  that  h(v)  be  corrected  to  the 
peak  height  that  a  Lorentzian  function  would  have  if  the  peak 
were  of  theoretical  width  (r  ).     This  eliminates  many  of  the 
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Table  1  —  Nomenclature  for  Experimental  Mossbauer  Spectrum  Parameters. 


Name  Symbol  Units 

Differential  do  cm/ s 

Chemical 

Shift 


Quadrupole 
Splitting 


cm/s 


Temperature  5^  cm/s/°K 

Coefficient  of 
the  differential 
chemical  shift 

Temperature  cm/s/°K 
coefficient  of  the 
quadrupole  splitting 

Line  width  F  cm/s 


Magnitude  e  percent 

of  the  Effect 


Definition 

Displacement  of  the  center  of  a  singlet  or 
a  multiplet,  from  the  center  of  the  doublet 
of  the  standard  reference  material. 

The  separation  of  two  peaks  resulting  from 
two  transitions  (1  and  2)  from  the  excited 
(prime)  to  the  ground  state  that  results 
from  electric  quadrupole  interaction  with 
an  electric  field  gradient,  where  m 
represents  the  magnetic  quantum  number 
resulting  from  this  interaction. 

The  variation  in  So  with  temperature  in 
degrees  Kelvin.  The  temperature  range 
should  be  specified. 


The  variation  of  AEq  with  temperature  in 
degrees  Kelvin.  The  temperature  range 
should  be  specified. 

The  full  width  at  half  maximum  of  the 
Mossbauer  peak. 


where  h('»)  =  height  of  resonance  base  line 
h(v)   =  height  of  peak  at  absorption 
maximum 
B   =  background 

h(v)  is  to  be  normalized  to  peak  height 

corresponding  to  theoretical  line  width. 


Asymmetry  S^  The  ratio  of  the  more  positive  velocity  peak 

Coefficient  height  to  the  more  negative  velocity  peak 

height  of  a  doublet.  For  normalization  see 

h(v)  above. 

Sg  The  ratio  of  the  more  positive  velocity 

absorption  area  to  the  more  negative  velocity 
absorption  area. 
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problems  of  peak  area  and  base  line  estimation.     It  is  neces- 
sary to  evaluate  the  background  which  results  from  such  non- 
resonant  radiation  as  Compton  scatter,  x-ray  f lourescence,  etc. 
Such  backgrounds  are  certain  to  be  different  for  different 
spectrometer  source-detector  configurations.     The  type  of  sam- 
ple matrix  will  also  alter  this  background  correction.     It  is 
to  be  understood  that  the  emphasis  is  on  direct  experimental 
parameters.     It  is  possible  to  relate  £  to  the  Debye-Waller 
factor  for  the  absorber  if  such  a  factor  is  also  known  for  the 
source . 

The  same  normalization  of  h  as  for  h(v)  described  above 

is  required  for  determining  S^.    The  asymmetry  coefficient 

which  measures  areas  (S^)  can  best  be  obtained  by  evaluating 

a 

the  function  parameters  of  the  Lorentzian  using  a  least  squares 
fit  to  the  experimental  data.    The  area  of  the  peak  can  then 
be  obtained  by  evaluating  the  function. 

An  example  of  a  spectrum  reporting  format  is  shown  on 
pages  kO  through  42.     The  first  page  contains  all  the  necessary 
chemical  data  such  as  purity.     References  on  synthesis ^  struc- 
ture^  and  other  work  on  the  spectroscopy  of  the  compound  are 
given.     On  the  second  page  is  the  spectrum  showing  a  best  least 
squares  fit  of  Lorentzian  functions  to  the  experimental  data 
points.     A  parabolic  baseline-  may  have  to  be  subtracted  as  well. 
For  transmission  spectra  the  peak  should  be  pointing  downward. 
On  the  third  page  are  the  calculated  parameters  and  their 
random  error  reported  as  a  standard  deviation  of  a  single  deter- 
mination.    All  of  the  unprocessed  data  in  the  spectrum  is  given 
on  this  page. 


39 


SODIUM  NITROPRUSSIDE 
DISODIUM  PENTACYANONITROSYL  FERRATE  DIHYDRATE 


Na   Fe(CN)  NO 

2  5 

SINGLE 


2H  0 

2 

CRYSTAL 


Supplier  —  Isomet  Corporation 


Muir 
Index  No. 

62C03 
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Phys.  Ctiem.  Solids 
23,  1783  (1962) 


Temperature 
Room 


8  cm/s 
-0.012   ±  0.005 


A  E  cm/s 
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L.  M.  Epstein,  J. 
Chem,  Phys.  36, 
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±  0.0008 
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±  0.0012 

W.  Neuwirth,  Z. 
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±  0.0003 

0.1725 
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Physik.  167,  176 

-25  °C 

-0.0592 

±  0.0003 
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25  °C 
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Phys.  41,  3378  (1964) 

P.  T.  Manoharan  and 
H.  B.  Gray,  J.  Am. 
Chem.  Soc,  87,  3340 
(1965) 

F.  A.  Cotton,  R.  R. 
Manochamp,  R.  J.  M. 
Henry,  R.  C.  Young, 
J.  Inorg.  Nucl.  Chem. 
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APPENDIX  D 


PHYSICAL  CONSTANTS  AND  CONVERSION 


FACTORS  FOR  IRON-57  MOSSBAUER  SPECTROSCOPY 


Nuclear  Parameters 


Excited  State  of 


57. 


Fe 


Energy,  E^  -  14.37  KeV,  0.864  A  [38] 
Half  life,  T^^2  =  9.8  x  lO^^s  [39] 


Internal  conversion  coefficient,  a  =  9.2  ±  0.5  [40,4l] 
Spin,  I  =  3/2"  [41] 

Magnetic  moment,  [i    =  -0.0154  nm  [39] 


Ground  State  of  -^'Fe 
Spin,  I  =  1/2"  [39] 

Magnetic  moment,  ^    =  +0.09024  t  0.00007  nm  [4l] 

Mossbauer  Parameters 

Cross  section,        =  I.9I  t  0.014  x  10~"^^cm^  [42] 

Natural  line  width,        =  4.655  x  10"^  eV  [39] 

=  0.009699  cm/s 


Quadrupole  moment,        =  +O.285  x  10' 


-24  2 


cm^  [393 


7.457  X  10"^-^  ergs. 


1.126  X  10^  hertz 
1.781  X  10"^^  kcal 
3.755  X  10"^  cm"-^ 


Iron  Atomic  Parameters 


Ground  state  configuration: 


x-ray  lines:  =  6.4  keV,  1.936  A  [43] 

=  7.06  keV,  1.75  A  [43] 
Mass  adsorption  coef f icient, -n/p  (l4„4keV)  = 

63.5  cmVg  [13,45] 
X-ray  scattering  cross  section,  a^(Fe^ l4 .4keV)  = 

4xl0"24cm2/atom  [44] 

Mossbauer  Iron-57  Standard  Reference  Material  Parameters 

o 

Quadrupole  splitting,  AE^  =  8.284  ±  0.0095  x  10~°  eV 

=  0.1726  t  0.0002  cm/s 
=  8.284  +  0.015  X  10"^^  ergs. 
=  2.003  ±  0.023  X  10'^  hertz 
=  1.979  t  0.036  X  10"^°  kcal 
=  6.682  +  0.077  X  10"'^  cm"-^ 

Percent  effect,   e  =  17^ 

Temperature  dependence  of  the  quadrupole  splitting, 

(Q^)  =  10-5cm,s-l,°  K"^[30] 

Temperature  dependence  of  the  chemical  shift,  6^^  = 

3  X  10-5cm,s-l,°K"^  [30] 

Gamma  attenuation  =  85^ 


For  -^'Fe 
1  cm/s  = 


only: 

4.798  X 
6.242  X 
4.136  X 
2.613  X 
1.240  X 


10"'  eV 
10^^  erg 
10-^5  hertz 
10^2  kcal 
10-4  cm"''' 
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APPENDIX  E 

ATTENUATION  AND  LINE  BROADENING  FOR  ^'^FE  ABSORBERS 


The  required  thickness  of  the  standard  reference  absorber, 
and  for  all  Mossbauer  absorbers.  Is  determined  by  the  gamma-ray 
attenuation  and  the  line  broadening  which  limits  the  resolution. 

The  attenuation  of  the  l4.4  keV  radiation  can  be  calculated 
from  the  mass  absorption  coefficient,  and  the  absorber  area 

density,  t^  (g/cm^),  using  the  equation: 

1  =  1  e"^^m^p 
o  ^ 

■while  for  compounds: 

U    t      =    E      LL  t 

^m  p      ^    ^m^  p^ 

The  subscript  referring  to  the  elements  forming  the  compound. 

The  mass  absorption  coefficients,  \i  ,  for  l4.4  keV  radla- 

m 

tion  plotted  as  a  function  of  element  [13]  are  shown  in  Figure  7. 

The  line  broadening  can  be  calculated  by  the  method  of 

Margulles  and  Ehrman  [33]  for  the  correction  of  Mossbauer  line 

width  for  source  and  absorber  thickness,  the  result  of  which  is 

shown  in  Figure  7.     The  width  can  be  found  from  the  effective 

absorber  thickness  T^  =  f '  t  ,  obtained  from  the  Debye-Waller 
a        a  p-^ 

factor,  f '     of  the  absorber,  and  the  absorber  thickness  in 
2 

g/cm    of  natural  abundance  iron.     Figure  8  is  a  plot  of  the 
broadening  of  the  transmitted  line  for  a  source  having  a  uni- 
form distribution  of  cobalt-57  atoms  on  a  non-ferrous  matrix. 
The  emission  and  absorption  spectra  are  assumed  to  be  Lorent- 
zlan.     It  is  to  be  understood  that  this  graph  is  to  be  used 
as  a  guide.     If  highly  reproducible  results  are  destired, 
mathematical  computation  should  be  used. 
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Figure  Attenuation  coefficients  for  l4.4  keV  gamma 

radiation  as  a  function  of  element. 
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EFFECTIVE  ABSORBER  THICKNESS 


Figure  8.     Plot  of  effective  source  thickness.  Adapted 
from  Margulies  and  Ehrman  [35]. 

f =  recollless  fraction  of  source  or  absorber 
where  subscript  l^a,  l=s  represents 
absorber  and  source  respectively, 
t     =  thickness  of  absorber  in  mg/cm^. 


P 


s     =  Mossbauer  cross  section  of  -^^ '  Fe  cm  /atom. 
Ng  =  number  of  ^Tpe  atoms  in  the  source  per  cm-^. 

t     =  thickness  of  source  in  cm. 
s 

a     =  experimental  full  width  at  half  maximum. 
=  effective  absorber  thickness  (mg/cm^). 

T    =  source  thickness. 
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Announcement  of  New  Publications  on 
Standard  Reference  Materials 

Superintendent  of  Documents 
Government  Printing  Office 
Washington,  D.  C.  20402 

Dear  Sir: 

Please  add  my  name  to  the  announcement  list  of  new 
publications  to  be  issued  in  the  series:    NBS  Miscellaneous 
Publication  260  -  Standard  Reference  Materials. 
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City  State   Zip  Code  

(Notification  key  N-330) 
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